Transition from band insulator to Mott insulator in one dimension: 
Critical behavior and phase diagram 
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We report a systematic study of the transition from band insulator (BI) to Mott insulator (MI) 
in a one-dimensional Hubbard model with an on-site Coulomb interaction U and an alternating 
periodic site potential V. We employ both the zero-temperature density matrix renormalization 
group (DMRG) method to determine the gap and critical behavior of the system and the finite- 
temperature transfer matrix renormalization group (TMRG) method to evaluate the thermodynamic 
properties. We find two critical points ai U — Uc and U = Us that separate the BI and MI phases 
for a given V . A charge-neutral spin-singlet exciton band develops in the BI phase {U < Uc) and 
drops below the band gap when U exceeds a special point Ue- The exciton gap closes at the first 
critical point Uc while the charge and spin gaps persist and coincide between Uc < U < Ua where 
the system is dimerized. Both the charge and spin gaps collapse at U — Us when the transition to 
the MI phase occurs. In the MI phase {U > Us) the charge gap increases almost linearly with U 
while the spin gap remains zero. These findings clarify earlier published results on the same model 
and offer new insights into several important issues regarding appropriate scaling analysis of DMRG 
data and a full physical picture for the delicate nature of the phase transitions driven by electron 
correlation. The present work provides a comprehensive understanding for the critical behavior and 
phase diagram for the transition from BI to MI in one-dimensional correlated electron systems with 
a periodic alternating site potential. 

PACS numbers; 71.30.-|-h, 71.10.Pm, 77.80.-e 



I. INTRODUCTION 



The nature of the insulating ground state of interacting 
electron systems has been a subject of long-standing in- 
terest and debate in condensed matter physics. Because 
of strong quantum effects caused by spatial confinement 
and technical advantages for theoretical treatment, one- 
dimensional (ID) electron systems have been most exten- 
sively studied^'^. Strong correlation effects in ID lead to 
the separation (decoupling) of charge and spin degrees of 
freedom. Starting from a gapless phase with charge-spin 
separated excitations, interactions can drive the system 
into new phases of different characteristics with (i) gap- 
ful charge excitations only, (ii) gapful spin excitations 
only, or (iii) co-existing gapful charge and spin excita- 
tions. For phases with both charge and spin excitations 
gapful, the charge and spin degrees of freedom are rarely 
decoupled. Despite these findings, there remain impor- 
tant unresolved issues regarding the quantum nature of 
the insulating state in ID interacting electron systems 
and the phase transitions driven by the electron corre- 
lation. The first issue concerns the establishment of an 
accurate phase diagram and the critical behavior near the 
phase boundaries. Secondly, a band insulator (BI) with 
quasi-particle excitations typically cannot be character- 
ized by charge and spin excitations. A proper characteri- 
zation scheme needs to be developed. Most importantly, 
the nature of the correlation-driven transition from BI to 
Mott insulator (MI) is still not fully understood. 



There has been considerable recent interest in the 
study of a prototype one-dimensional model for ferro- 
electric perovskites for the understanding of the response 
of strongly correlated electron systems with lattice dis- 
tortions. These efforts have raised and addressed some 
fundamental issues in the nature of the quantum phase 
transition and the related critical behavior in ID interact- 
ing electron systems. Earlier works^iii^ mainly deal with 
the effects of strong electron correlation on the electron- 
lattice interaction and the polarization effects in the in- 
sulator. Quantum phase transitions and the character- 
ization of the insulating state are the focus of more re- 
cent work.^i'^i^i^i^°i"i^'^i"i^^i^^i^^ In particular, an issue 
of fascinating debate is the nature of the transition (or 
crossover) from the BI phase to the MI phases. Gidopou- 
los et al. showed^^ that due to the reversal of inversion 
symmetry of the ground state from BI to MI, there is a 
critical point for spin excitations. However, for charge 
excitations the critical behavior is less clear. Recently, 
Fabrizio et ali^ developed an effective field theory for 
this problem and showed that there are two continuous 
transitions from BI to MI . One is a spin transition of 
the Kosterlitz-Thouless type at a critical point U — Us, 
and the other is an exciton transition at an Ising critical 
point U = Uc < Us where the exciton gap closes. Be- 
tween Uc and Us , the site-parity is spontaneously broken 
and the system is characterized by a doubly degenerate, 
dimerized ground state. These results raise interesting 
questions about the structure of the ground-state phase 
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diagram of ID interacting electron systems and the char- 
acterization of the critical behavior near the transition 
points from the BI phase to the MI phase. 

The model Hamiltonian for the system of interest is 
defined in the Hubbard formalism at half-fillingS*Si2ii£ 



Ac(L), and (iii) the spin-triplet gap As(L), 



H = 



J2 [-t {4.c^+la + h.c.) + V{-iyn, 

i<7 



2 ' 



(1) 



where and are the electron creation and number 
operators at site i, U > is the on-site Coulomb repul- 
sion, and V is the staggered site chemical potential. This 
model captures the key ingredients of one-dimensional 
correlated insulators with mixed ionic-covalent charac- 
ters, such as oxide dielectric materials^" and quasi-lD 
organic charge-transfer complexes. It incorporates co- 
valency, ionicity, and strong electron correlational^ 

In this paper, we present the results of exten- 
sive calculations for Hamiltonian (1) using both the 
zero-temperature density matrix renormalization group 
(DMRG) methocS and the finite-temperature transfer 
matrix renormalization group (TMRG) method^^i^^iS^. 
These methods have been demonstrated to be highly ac- 
curate for ID interacting electron systems. Our aim is to 
systematically examine and clarify issues raised in recent 
work and to provide a comprehensive understanding for 
the transition from BI to MI in one dimension. We show 
detailed results on the gap and critical behavior. We find 
that a charge-neutral spin-singlet exciton band forms in 
the BI phase and drops below the band gap as U in- 
creases beyond a special point Ue- With increasing U the 
excitons then condense and the system enters a dimer- 
ized phase, followed by the closure of the quasi-particle 
(both spin and charge) gap when the system enters the 
MI phase. We clarify basic concepts on charge and spin 
excitations studied in recent work. Our results support 
the conclusion of Fabrizio et al}^ on the existence of two 
critical points for the transition from BI to MI. We also 
present detailed results on the formation of the exciton 
band and the scaling behavior near the critical points. 
In addition, we carry out TMRG calculations to study 
thermodynamic properties to further elucidate the gap 
and critical behavior of the system. 



II. LOW-ENERGY EXCITATIONS 

To properly characterize the BI and MI phases and es- 
tablish the phase diagram, we need to evaluate the behav- 
ior of several low-energy excitations, including the charge 
and spin excitations and an exciton excitation that will 
be used to characterize the BI phase. We calculate the 
following three excitation gaps defined on a finite ID lat- 
tice of length L (chosen as an even integer) at half-filling: 
(i) the singlet exciton gap Ae(L), (ii) the charge gap 



Ae(L) 
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where {N^ ,Ni) is the lowest energy of the system with 
7V| up and down spin electrons, and Ei{Ni, N^) is the 
lowest energy of the singlet excitations. For the charge 
gap Ac, there is an alternative definition 



Ae(i) = Ei 



1, 



L 



1) 



r2' 



(5) 



For the model studied in this work, these two definitions 
on the charge gap give the same results in the thermody- 
namic limit. This is supported by our numerical calcu- 
lations. However, since the numerical accuracy is much 
higher in calculating Ac by using Eq. (PJ than by using 
Eq. 0, wc use Eq. Q in all the reported calculations. 

Although Ac and Ag are usually considered to be the 
charge gap and the spin gap, respectively, in the litera- 
ture, Ac in fact is the chemical potential jump for par- 
ticles in the system. It measures the chemical potential 
jump of putting a particle into or taking a particle out of 
the system. Only when the charge and spin excitations 
are separated and the spin gap is zero, Ac is equal to 
the charge excitation gap as, for example, in the stan- 
dard Hubbard model (Eq. (P) with = or [/ — > oo). 
When the charge and spin excitations are not separated, 
or when the spin gap is not zero, the chemical potential 
jump is not equal to the charge excitation gap. 

At half-filling, the first excitation state can be either 
a spin singlet or a spin triplet and Ac never exceeds 
Ag. When Ag < A^, the first excitation state must be 
a charge-neutral spin-singlet state; otherwise, Ac equals 
As and measures the excitation gap of an exciton band. 

In the presence of a nonzero V in the Hamiltonian, 
the MI phase is reached when U is large (including the 
limit U oo). In the MI phase, it is well understood 
that the charge gap is non-zero, but the spin gap is zero. 
Therefore Ac is the charge excitation gap. Meanwhile, 
both As and Ac are zero and the gapless elementary 
excitations in the system are spinonsi^ 

In the BI phase, all the elementary excitations are gap- 
ful. Let uj.^ and dj.^ be the creation operators for parti- 
cles in the upper conducting band and for the holes in the 
lower valence band. The Hamiltonian for free particles 
at J7 = is 



H 



k(7 



(6) 



where k is the momentum, a is the spin index, and 
Ek = VV^ + 4i2 cos2 k. We have Uka\GS) = vljGS) = 
for the ground state \GS) at half filhng with Uka be- 
ing the annihilation operator for electrons in the conduc- 
tion band, and the creation operator for electrons 
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in the valence band. For particle-hole excitations with 
one particle and one hole in the system, it is clear that 

When U is small in the BI phase, the particle or hole 
excitations become dressed quasi-particles, and Ac mea- 
sures the chemical potential jump of the particles. Since 
the charge and spin degrees of freedom are not sepa- 
rated, Ac is not exactly the "charge gap" derived from 
the gapless charge-spin separated excitations in the field 
theoretical approach. If the particles and holes are not 
bounded, then Ac = A^ = Ag. However, the particle- 
hole excitations may bound to form excitons, and result 
in an exciton gap smaller than Ac. It is confirmed by 
our DMRG calculations shown below that Ae equals Ac 
when U is small, but becomes smaller than Ac when U 
exceeds a special point C/e, indicating the formation of 
singlet excitons. Meanwhile, the calculations show that 
As = Ac is always observed, indicating that there is no 
triplet exciton formation in the BI phase. 

The DMRG calculations also show that in the BI phase 
the exciton gap and quasi-particle excitation gap de- 
crease and approach zero with increasing U . The ex- 
citon gap closes at the first critical point Uc, and the sys- 
tem enters a phase where the excitons condense to form 
dimerized ground state. At the second critical point Us 
the quasi-particle gap also closes (Ag = A^ = Ac = 0). 
When U exceeds Us, Ag and A^ remain zero but Ac in- 
creases almost linearly with U . Therefore, Us is the point 
where the quasi-particle excitation gap in the BI phase 
collapses, and spinous in the MI phase form. This pic- 
ture is consistent with that of the recent field theoretical 
studiesi^. 

From the physical picture outlined above, it is clear 
that there are three special points Ue, Uc, and Us along 
the U scale. Among them Uc and Us are two critical 
points separating the BI and MI phases while Ue is a 
special point signaling the formation of the spin-singlet 
exciton band in the BI phase. In the thermodynamic 
limit for a given V the system can be divided into the 
following regions: 

1. < U < Ue, Ae = Ac = As ^ 0, gapful quasi- 
particle excitations only. 

2. Ue < U < Uc, < Ac < Ac ^ As, gapful quasi- 
particle excitations coexist with singlet particle- 
hole excitations that bound to form excitons. 

3. U = Uc, Ae = 0, Ac = As > 0, the system is 
critical for exciton excitations. 

4. Uc < U < Us, Ae < Ac — As, excitons condense, 
and the system is dimerized. 

5. L/ = t/s, Ac = Ac = As = 0, the system is critical 
for quasi-particle excitations. 

6. U > Us, Ae — As ^ 0, Ac > 0, gapless spinon ex- 
citations in the MI phase; the particle-hole picture 
breaks down. 



In the following we present the calculated results on 
the gap and critical behavior leading to the establishment 
of the phase diagram. In all reported calculations, free 
boundary conditions are used in the zero-temperature 
DMRG calculation. For a finite size system, we can show 
rigorously using the variational principle that the lowest 
energy state of the Hamiltonian with free boundary 
conditions in each (N-^jNi) subspace is non-degenerate 
except for an up-down spin degeneracjtSS. Therefore 
there is no level crossing with the lowest energy state 
in the {Ni,Ni) subspace and i?o(A^t, A^j.) is an analytic 
function of U and V. From this property and Eqs. Q 
and Q), we can further show that Ac(-L) and As{L) are 
also analytic functions of U and V. In this work we fo- 
cus on two values of the staggered potential V=0.3t and 
l.Ot and study the behavior of the three gaps introduced 
above in response to the on-site Coulomb repulsion U. 



III. EXCITATION GAPS 

For fermion systems, the truncation error of DMRG 
iterations is generally much smaller than that of spin 
systems when the same number of optimal states are re- 
tained. The efficiency of the finite system DMRG method 
is related to the truncation error; the bigger the trunca- 
tion error is, the larger the improvement of the finite 
lattice sweeping can make. We used both finite and in- 
finite lattice DMRG algorithms in testing calculations. 
We find that the improvement of the ground-state en- 
ergy made with the finite lattice sweeping is very small 
when a large number of states are retained. A better way 
to increase the accuracy of the results is using the infinity 
lattice approach by retaining more states. 

Figure n shows the behaviors of As, Ac and Ac as a 
function of 1/L for different t7 at ^ = 1.0^. For V = OM, 
similar results can be drawn but the V — l.Ot case is more 
accurate because most of the features can be seen when 
the chain length is short, while in the V — 0.3t case, 
very long chains need to be used to obtain the same re- 
sults. Figuren(a) presents the results for U = 2.0t and 
clearly shows that the three gaps converge to the same 
finite value in the thermodynamic limit. The difference 
between these gaps shown in Fig. ^ (b) also displays this 
feature clearly. When the chain length is short, the ex- 
citon gap is larger than the spin gap, and level crossing 
happens at a finite chain length where the exciton gap 
drops lower thereafter. The exciton gap decreases con- 
tinually and reaches a minimum when the chain length 
increases further; it then starts to increase and converge 
to the value of the spin and charge gaps. For U — 3.0t, in 
Fig. n (c) and (d) , the spin gap and charge gap still con- 
verge to the same value in the thermodynamic limit, but 
the exciton gap goes to a different value lower than the 
other two gaps. In short chains, the exciton gap is still 
larger than the spin gap, but after they cross each other, 
the exciton gap decreases monotonically. One can also 
see that the second state in the singlet sector also crosses 
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FIG. 1: Energy gaps vs 1/L for V = l.Ot and (a) V = 2.0t, 

(c) U = 3.0t, (e) [/ = 4.0t. Solid lines are Ae, dot-dashed 
lines are As, and short-dashed lines are Ac. The long-dashed 
line in (c) represents the gap of the second state in the sin- 
glet sector Ac^'. The dotted lines in (a) and (c) denote the 
thermodynamic limit of the corresponding gaps. Panels (b), 

(d) , and (f) show the absolute values of the difference of the 
charge and exciton gap with the spin gap, where empty circles 
show Ac — As and filled circles show Ae — As. 500 states are 
retained in the DMRG calculations. 



the spin and. charge gap and converges to the lowest state 
when L — > 00. In our calculations, we also see that more 
states in the singlet sector cross the spin and charge gap 
with increasing chain length. This shows that the whole 
spectrum of the exciton sector decreases in value and the 
exciton gap is indeed different from the other two gaps. 
For the case oiU = 4.0t, shown in Fig. ^ (e) and (f), 
no level crossing for different chain length is detected. 
All three gaps decrease monotonically when the chain 
length increases. The exciton gap and the spin gap ap- 
proach zero at infinite chain length, indicating that there 
is no gap for the exciton and spin sectors. Meanwhile 
the charge gap approaches a finite value in the thermo- 
dynamic limit. 

For all the cases we have studied, Ac decreases mono- 
tonically with increasing L. However, the size depen- 
dence of As is more complicated. In certain ranges of 
U jt and Vjt close to the critical regimes, including the 
case shown in Fig. ^ (a) and (c) , A^ and Ae vary non- 
monotonically and their minima are located at a finite 
L = imin rather than aX L — 00. In a recent work^S^ 
the authors studied the same model Hamiltonian ^ us- 
ing the DMRG method but did not observe such non- 
monotonically behavior, and suggested that such a be- 
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FIG. 2: The spin gap As vs 1/L for V = IM and U = 2.50t. 
The number of retained states are m=300 (empty circles), 400 
(empty squares), 500 (filled circles), and 800 (filled squares). 



havior may be due to the loss of accuracy in DMRG 
calculations when the chain length is increased or due 
to some intrinsic length scale for the spin degree of free- 
dom. We have carefully examined this issue by carrying 
out extensive scaling analysis. We demonstrate that the 
non-monotonical behavior is not due to the lack of ac- 
curacy of the calculations, instead the behavior is a true 
feature of Hamiltonian ^ with the open boundary con- 
dition (OBC). Figure 121 shows the chain length depen- 
dence of the spin gap for U — 2.5t and V — l.Qt cal- 
culated by retaining different numbers of optimal states 
m = 300,400,500, and 800. One can see that the min- 
imum occurs at L ^ 30; at this length the accuracy of 
the DMRG calculations are still very high. More signifi- 
cantly, the results for different m fall onto the same curve 
(except for the cases of L > 100 and m = 300). It shows 
unambiguously the existence of the minimum of As in 
its dependence on the chain length. For the exciton gap 
Ae, the situation is the same. In fact the occurrence of a 
gap minimum at a finite L is not an uncommon feature 
for a system with incommensurate low- lying excitations. 
It suggests that the spin excitations of the model Hamil- 
tonian maybe incommensurate with a characteristic 
wave vector defined by 27r/Lmin (or tt — 27r/L,„in) in some 
area of the phase space. 

Comparing Fig. ^ (a) and (c) , it is clear that there is 
a special point Ue, where the exciton gap begins to devi- 
ate from the spin gap (for V = l.Ot, 2.0t < Ue < 3.0t). 
In both cases, all three excitations are gapful. The sys- 
tem is in the same (BI) phase as the U = case, where 
As = Ac = Ae = 2V in the thermodynamic limit. 
Here particle-hole excitations bound into excitons by the 
Coulomb interaction at C/ > Ug. Fig. O shows the dif- 
ference As - Ae for V = l.Ot and 0.3t. The fitting to 
DMRG results gives the critical value Ue — 2.26At for 
V = l.Ot and Ue = 1.276i for V = 0.3t. 

For finite L we find that Ac is always larger than As. In 
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FIG. 3: The difference between tfie spin gap and tfie exciton 
gap near Ue- (a) V=f .Ot, 2.264t. Tlie daslied fitting line 

is 0.026([//f - 0.ni){U/t - 2.264); (b) V=0.3t, (/, - f .276t. 
Tlie dasfied fitting line is Qm9>2[U/t - 0.584)((7/f - f.276). 



the MI phase, Ac is finite but approaches zero in the 
thermodynamic limit. In the BI phase, A^ and Ac always 
approach the same value in the thermodynamic limit. 
This can be seen either from the asymptotic behaviors of 
As and Ac in the limit L — > oo [Fig. ^{a) and (c)] or from 
the 1/L dependence of the difference Ac — As (Fig. ^(b) 
and (d)). For all the cases we have studied, we find that 
Ac — As drops monotonically and approaches zero in the 
limit 1/L — > even when As changes non- monotonically. 
For a given V, this result holds from f/ = up to a 
critical regime where both Ac and As become smaller 
than truncation errors. It suggests that Ac and As are 
equal in the thermodynamic limit in the entire BI phase. 

When As changes non- monotonically with L, the ex- 
trapolation for the spin gap in the limit 1/ L becomes 
subtle. If the data with L < L^in are used in the extrap- 
olation, the extrapolated value of As will certainly be 
smaller than the true value. However, if the data with 
L > Lmin is used but L is still not large enough to reach 
the regime where As begins to saturate, the extrapolated 
value of As will be larger than the true value (this seems 
to be the case in Ref. |23). For the data shown in Fig. 
^a) and (c), these two kinds of extrapolations result in 
Ac > As and Ac < As, respectively. Both are incorrect. 



To correctly extrapolate As in the limit 1/L ^ 0, data 
with L much larger than Lmin must be used. 

For the U — A.Ot case shown in Fig. ^(s), Ac is finite 
but As becomes zero in the thermodynamic limit. This 
indicates that the system is in the same phase as that 
for U oo, namely the MI phase. In the MI phase, the 
chain length dependence of the three gaps are monoton- 
ical; in addition, there is no crossing between As and Ac 
when the chain length varies. In the L — > cx) limit the 
spin gap As is zero in the MI phase suggesting that there 
is a critical point that separates the MI phase from the 
BI phase. At this critical point Us, the spin gap vanishes. 
Because the charge gap is equal to the spin gap in the BI 
phase, the charge gap will also vanish at the same point 
Us- However, when U increases further, the charge gap 
increases with U while the spin gap remains zero. Con- 
sidering that the exciton gap is lower than the spin gap 
in the BI phase at [/ > C/g, Ag may vanish before the 
spin gap and charge gap do. In that case, there should 
be another critical point Uc signaling the collapse of the 
exciton gap. 
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FIG. 4: The spin (empty circles), charge (solid circles) and 
exciton (empty squares) gaps as a function of U for V = 
O.St (a) and V = l.Ot (c). Panels (b) and (d) show the U 
dependence of the spin and exciton gaps in the vicinity of the 
critical region for (a) and (c) respectively. At Ue, the exciton 
gap deviates from the other two gaps and it collapses first 
at Uc - 2.225t for V = 0.3t and ~ 3.675* for V = l.Ot. The 
charge and spin gaps collapse at Us ~ 2.265t for V — 0.3t and 
~ 3.71t for V = l.Ot. For < U < Us the spin and charge 
gaps have the same value in the thermodynamic limit. When 
U > Us, the spin and exciton gap are zero while the charge 
gap is finite. 

In Fig. 21 we show the U dependence of the three gaps 
for V = 0.3t [(a) and (b)] and V = l.Ot [(c) and (d)]. For 
both cases, there are indeed two critical points Uc and 
Us although they are very close. When U < Ue, Ae — 
As = Ac, and the three gaps decrease almost linearly 
with increasing U. At Ue, the exciton gap splits off and 
drops below the other two gaps. At the critical point Uc, 
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the exciton gap collapses while the spin and charge gaps 
still coincide and remain finite until the second critical 
point Us where they both collapse. At J7 > Us, the charge 
gap increases with increasing U while the spin gap and 
the exciton gap remain zero in the thermodynamic limit. 
Although the accuracy of our DMRG calculations do not 
allow a direct assessment of the behavior of the exciton 
gap for Uc < U < Us, we believe that Ae is finite in this 
region (see more detailed discussion on this point in the 
following section) . The extrapolation of the gap behavior 
leads to U^ ^ 2.22U and Us ^ 2.265t for V = 0.3t, while 
Uc 3.675t, Us - 3.71t for V = IM. 

It is clear that the U dependence of the three gaps is 
similar for V — 0.3t and V = l.Ot. It is expected that the 
same picture is valid for all V. Furthermore, Ue, Uc, and 
Us all approach the same point Uoo = "^V in the ^ oo 
limit. 



IV. CRITICAL BEHAVIOR 
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FIG. 5: The behavior of the charge gap in the vicinity of 
17s for V = 0.3t. (a) The chain length dependence of the 
charge gap at [/ = 2.2t for different numbers of optimal states 
retained: m—200 (filled circles), 250 (empty squares), 300 
(empty circles) and 400 (filled squares) ; (b) The dependence of 
the charge gap on the number of the retained states at L=300 
with U = 2.2t; (c) The charge gap for m — > oo with different 
chain lengths: L — 300 (empty circles), 400 (filled circles), 500 
(empty squares) and 600 (filled circles) in the vicinity of the 
critical region; the fitting lines using eq. (|HJ are also shown; 
(d) The chain length dependence of Ac,min{L); (e) The chain 
length dependence of Ua,L/t; (f) The chain length dependence 

of Oc.L. 

From the gap behavior presented in the previous sec- 



tion, it is clear that there are two critical points Uc and 
Us for Hamiltonian for a given V. It is important 
to study the detailed critical behavior near the critical 
points for the understanding of the nature of the BI- 
to-MI transition. In the following we study the critical 
behavior of the system by examining (i) the evolution 
of the gap behavior near the critical points and (ii) the 
behavior of the ground-state energy of the system. 

A. Analysis of the gap behavior 

The U dependence of the gaps shows that the charge 
instability occurs at the same point as that for the spin 
transition. To determine the critical points for the charge 
and spin excitations, we examine when Ac and A^ be- 
come zero in the limit L ~* oo. Since the numerical errors 
are larger than the magnitude of Ac or As in the vicinity 
of the critical points, it is difficult to determine accu- 
rately the critical behavior simply from the values of the 
energy gaps. To resolve this issue, we analyze the scal- 
ing behavior of Ac{L) around its minimum with respect 
to U. However, the DMRG results of Ac{L) depend the 
number of states m retained during the iterations. When 
the chain length is long enough, the difference due to re- 
taining different number of states show clearly. In Fig. [S] 
(a), we show the chain length dependence of the charge 
gap at U = 2.22/: and V = 0.3t by keeping different m. 
The difference is obvious. This problem can be solved by 
employing the extrapolation in the limit of m — > cxd. For 
given U and V, the charge gap at chain length L and by 
keeping m states is Ac(to, L). By extrapolating to the 
infinite m limit, more accurate result of the charge gap 
at chain length L can be obtained, 

Ac(oo,L)= lim Ac(m, L). (7) 

m — *oo 

In Fig. O (b), we show the 1/m dependence and the 
extrapolation procedure of the Ac(rn, L) for L=300 at 
U = 2.22t and V = 0.3t. The extrapolated resuh 
Ac(oo,L) is considered the exact charge gap Ac{L) at 
chain length L. 

The extrapolated charge gap Ac(oo, L) is a function of 
U and chain length L at a given V. By applying the same 
procedure shown in Fig. 13 (b), we obtain Ac(oo,L) for 
different values of U near the critical point Us for a serial 
of selected L. In Fig. EKc), we show the U dependence of 
Ac{L) (Ac(oo,i)) in the vicinity of the critical point Us 
at chain length L — 300, 400, 500, and 600, respectively, 
for V^=0.3t. A gap minimum at finite chain length is 
clearly seen. Assuming Ac,min(i) to be the minimum 
of Ac{L) located at Us.l, then around this minimum we 
can expand Ac{L) to the leading order of the parameter 
u = U — Us,L as 

Ac{L) = Ac,,„i„(L) -t- ac {L) + O (u^) . (8) 

Since Ac{L) is an analytic function of U, both Ac,min(-^) 
and ac (L) should be finite. The critical behavior of 
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the charge excitations is determined by the properties of 
^c,min(i) and Uc (L) in the hmit L — > oo. If Ac_min(-^) 
in the hmit L ^ oo, the charge excitation is critical at 
Us = Us,oo, which would be consistent with the discus- 
sion in the previous section. However, if Ac,min remains 
finite in the limit L —t oo, then there is no critical point 
for charge excitations and the ground state is insulating 
in the entire parameter space. 

Figure[Sl(d) shows the calculated Ac,min(^) as a func- 
tion of 1/L. The solid curve is a least-square fit of the 
data and given by Ac,min(i) ~ 17.894/L - 729.785/L2. 
Within numerical errors, we find that Ac,min(^) is in- 
deed in the limit 1/L — > 0. Figure (e) shows the 
L dependence oi Us,l- It changes almost linearly with 
1/L. Within numerical errors, we find that the data of 
Usx are well fitted by Us,L/t « 2.265 - 22.532/L. Thus 
the critical point is at Uc/t = 2.265t, in full agreement 
with the value obtained by fitting the gap directly in the 
previous section. 

We now turn to the critical behavior of ac [L). Fig- 
ure (f) shows the 1/L dependence of ac{L) for the 
case V — Q.it. The fitting curve (solid line) is given by 
— 0.323-1-0. 0091L. The divergence of ac (L) suggests that 
the derivative of Ac{L) is singular at Us and the leading 
term in Ac in the thermodynamic limit is linear rather 
than quadratical in u, i.e., Ac{U) ^ \U — Us\- 



B. Analysis of the ground-state energy 

The ground-state energy as the zero-temperature free 
energy can also provide evidence for the critical behav- 
ior. However, singularities in the ground-state energy 
are of higher order derivatives with respect to the model 
parameter for continuous phase transitions. As a result, 
evidence for critical behavior derived from the ground- 
state energy is not as strong as that from the gap behav- 
ior, this despite the higher accuracy of the ground-state 
energy than that for the gap in the DMRG calculations. 

We have calculated the ground-state energy in the crit- 
ical region by retaining m=800 states and up to chain 
length L=1000. For Hamiltonian with open bound- 
ary conditions, the ground-state energy per site eo{L) 
satisfies 



eo{L) = 



Eo{L ) 
L 



= '^ + Y + j2+o{j^)^ (9) 



here Eq{L) is the ground-state energy for a chain of 
length L, eo the ground-state energy per site for L oo, 
Cb is the boundary energy (surface energy) due to the 
free boundary condition, and c — vtt where v is the spin 
wave velocity. When U > Uc, the system is gapful and 
c should approach zero when the chain length is much 
larger than the correlation length ^. 

The V = 0.3t, L > 200 ground-state energy results are 
fitted directly by 



00 









_ °°o (b) 

O 




O 

\ 


rP 





(d) 



40 
20 

< 
-20 
-40 



-1.2 
15 



pq 



2.24 
U 



2.24 
U 



eo{L) = eo + Cb/L + c/L"^, 



(10) 



FIG. 6: (a), (b), (c): The U dependence of eo, and c 
obtained from fitting tlie L > 200 data of the ground-state 
energy; (d) The chain length dependence of c(L) for U = 1.5t 
(dashed line), U = 2.0t (dot-dashed line), and U — 2.2t (solid 
line); (e) A = 10®(eo(L) - eo(lOO)) for L=120, 150, 200, 250, 
300, and 400 (from bottom to top on the ^ > side); (f) 
B — 10^(e(,(I/) — 66 (100)) for the same chain lengths as in (e) 
(from top to bottom on the B <Q side). 



and the obtained results are shown in Fig. El(a), (b) and 
(c). eo and Cf, are analytic functions of U . For c. Fig. (B^c) 
shows that it is not only non-zero for U < Uc but also 
has a fairly large value. This is due to the finite chain 
length effect. The value of c depends very sensitively on 
the chain length range used for fitting. 

To analyze the chain length dependence of c, we fit the 
ground-state energy for L — 20, L, and L -I- 20 using Eq. 
((Tn|l and vary L from 24 to 980. Each fitting gives the 
exact result of eo{L), eb{L) and c{L). In Fig. Efd), we 
show the L dependence of c(L) for U = 1.5t, U — 2.0t 
and U = 2.2t. It is clear that for U = 1.5t and 2.0t, c{L) 
becomes zero when L ^ oo. The result oiU = 2.Qt shows 
a minimum at L ^ 90, and the c(L) begins to approach 
zero after the minimum. For U = 2.2t, a minimum is 
also clearly seen. A comparison of these results leads 
to the conclusion that for L oo, c{L) will approach 
zero. The largest chain length used in the calculation is 
not long enough to obtain correct c(L) values. However, 
the obtained £o{L) and eb{L) may display the critical 
behavior of Uc and Us. In Fig. El(e) and (f), we show the 
results of e.o{L) — eo(io) and eb{L) — eb{Lo) for Lq = 100. 
These results show the existence of the critical point Uc. 
Lq can be viewed as a characteristic length of the critical 
region. 

From the U dependence of eq, it is possible to examine 
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the type of transition at the critical points. A problem 
is that the fitting results shown in Fig. El^a) include ex- 
tra errors induced by the fitting method. To avoid this, 
we analyze the U dependence of the ground-state energy 
using a different approach. 
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FIG. 7: (a) d^eo/dU^ versus [/ for L = 500 (filled circles), 600 
(filled squares), 800 (empty circles), 1000 (empty squares). 
The solid line is an extrapolation of the data to the limit 
L —> oo. Inset: Enlarged figure in the vicinity of the critical 
point Us- (b) The dependence of the d^eo/dU^ with ln(l/I/) 
at U = 2.21i (empty squares), 2.225t (filled circles), 2.27t 
(empty circles). The solid fitting -0.00997+0. 0091n(l/L) is 
the least square fitting for the U = 2.225t case. 



five does not diverge. These observations suggest that 
the phase transition at Uc is of the second order. No 
singularity is found in the first derivative or the second 
derivative near the critical point Us- This means that 
the transition at Us is higher than second order. These 
results are consistent with those reported by Fabrizio et 



V. TMRG STUDY OF SPIN SUSCEPTIBILITY 
AND SPECIFIC HEAT 

To gain more insight into the physics of the BI-to-MI 
transition, we have studied the thermodynamic proper- 
ties of the model using the TMRG method^^i^'^i^^ which is 
implemented in the thermodynamic limit and can evalu- 
ate very accurately the thermodynamic quantities at low 
temperature for quasi-lD systems. In our calculations, 
we kept 250 optimal states. The calculated specific heat 
C^, charge susceptibility Xc, and spin susceptibility Xs 
for U/t = 1.0,2.25,5.0 and V/t = 0.3 as a function of 
temperature are shown in Fig. |S1 




At a finite chain length L, the ground-state energy per 
site eo(i) is also a function of U. Here eo{L) contains 
only the errors from the DMRG truncation. We can ex- 
amine the derivatives of eo{L) with respect to U and 
analyze their chain length dependence. In Fig. Cfa), we 
show the second derivative of eo{L) with U for L=500, 
600, 800, and 1000. At each chain length, there is a min- 
imum near Uc in the U dependence of the second deriva- 
tive. When the chain length increases , the position of 
the minimum moves towards larger U and approaches Uc 
which is the critical point in the thermodynamic limit; 
meanwhile, the shape of the minimum becomes sharper. 
Fig. I3b) shows the chain length dependence of the sec- 
ond derivative. It is clear that for U = 2.225t ^ Uc, 
the second derivative diverges logarithmically with the 
chain length, but for other values of U the second deriva- 
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FIG. 8: Temperature dependences of (a) the spin susceptibil- 
ity Xs, (b) the charge susceptibility Xc, and (c) the specific 
heat Cv (t is set to 1). 

We find that Xc decreases exponentially at low tem- 
peratures in both the BI and MI phases, while Xs shows 
activated behavior only in the BI phase. In the MI phase, 
there are two broad peaks in C^, probably due to the 
charge-spin separation. Near the critical point, U=2.25t, 
since both the charge and spin energy gaps are very small, 
the exponential decays in both Xc and Xs show up only 
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at very low temperatures. These results support the con- 
clusions of the DMRG calculations presented in previous 
sections. 



VI. PHASE DIAGRAM 

The overall U dependence of the charge, spin and exci- 
ton gaps shown in Fig. 0]give a lot of information on the 
phase diagram of Hamiltonian . The charge and spin 
gaps coincide in the BI phase. Above Us, Ac increases 
with U but As remains zero. The exciton gap Ag col- 
lapses at Uc < Us, and when U > Us, the exciton gap 
should also be zero. However, from Fig. ^it is unclear 
whether the exciton excitations are gapful or gapless in 
the regime between the two critical points, Uc < U < Us- 
Even if the exciton excitations are gapful in this regime, 
the gap would be too small to detect numerically. 
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Here we evaluate the dimer- 



2? = ^E(-l)'(4cm. + /i.c.). 



(11) 



Figure Ela) shows the chain length dependence of the 
dimerization operator for different U. It is clear that 
for U = 2.2 li, when L — > oo, I? approaches zero. At 
U = 2.29t, V just starts to fall at the largest chain length 
we studied; it is expected that it will approach zero as 
the chain length is long enough. For U — 2.25t, which 
is between the two critical points, it seems that V will 
diverge to a nonzero constant. For a large range of chain 
lengths, the results can be well fitted by a straight line 
shown in Fig. 1^1 (a). The dependence of the finite chain 
dimerization on U is shown in Fig. Elb) for L=200, 
300, 400, and 500. These results indicate that in the 
thermodynamic limit, the ground states are dimerized 
when Uc <U < Us- 

The dimerization of the ground state for Uc < U < Ug 
suggests that the exciton excitations are gapful in this 
region. So the physical picture on the exciton excita- 
tion is emerging: the exciton gap formed in the BI phase 
collapses at the critical point Uc', with further increas- 
ing U, the (small) exciton gap will first increase, reach a 
maximum and then decrease and collapse again at Us', at 
U > Us, the exciton excitations remain critical. 



2 - 








FIG. 9: (a) Chain length dependence of the dimerization or- 
der parameter T> for U = 2.21t (empty squares), U — 2.25t 
(filled circles) and U = 2.29t (empty circles). The straight fit- 
ting line is 0.0499 + 3.886/L. (b) I){L) in the critical regime 
for L=200 (empty circles), 300 (filled circles), 400 (empty 
squares), and 500 (filled squares). The dotted lines indicate 
the two critical points. 

When the exciton gap collapses, the excitons can con- 
dense into the ground stataSS,. In this case, the system is 



FIG. 10: The ground-state phase diagram for Hamiltonian 
Q. The empty circles denote the DMRG results. The curves 
of Us and Uc are shown with solid lines that are very close to 
each other. The curve of Ue is shown by the dashed line which 
does not indicate a phase transition line, but rather denotes 
a serial of special points. The dotted line is U = 2V which is 
the limit for V — > oo. 

When U <^V, first order perturbation leads to 

As{U) = Ac{U) ^ V - cU, (12) 
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where c = V d k{y/2TTek)~^ is a constant deter- 
mined by the single-particle energy dispersion Sk = 
y/V^ + cos^ k. Since Ac{U) drops almost linearly 
with U/t in the BI phase, we can estimate the value of 
Us from Eq. l(T^ as Us ~ V/c. In the limit F — » 0, we 
have 



V/t « cie-^^*/^= 



or 



Us/t: 



iog(y/cit) ' 



(13) 



where C\ and C2 are two constants of order one. Figure lTUI 
shows the ground-state phase diagram for Hamiltonian 
(Q. The curve for Us and V < l.Ot is obtained from 
Eq. (|13|l . The parameters Ci and C2 are fixed by the two 
Us values for V = OM and V = l.Ot. When V 0, 
Us goes to zero but the ratio Us/V diverges. In the limit 
U/t ^ oo, Us is very close to 2V. The difference between 
Us and 2V is of order t: Us ~ 2V t. 



eral important issues raised in recent works and have es- 
tablished the ground-state phase diagram. We have iden- 
tified two critical points, Uc and Us, that separate the BI 
and MI phases. When U > Us, the system is in the MI 
phase where the charge excitations are massive but the 
spin excitations are critical. When U < Uc, the system 
behaves like a classic band insulator: the charge and spin 
excitation gaps coincide and a charge-neutral spin-singlet 
exciton band forms below the band gap when U exceeds 
a special point Ug. Between the two critical points, exci- 
tons condense and the ground state is dimerized. These 
results are consistent with the conclusions of a recent field 
theoretical study of the same model. The present work 
provides a detailed account for the critical behavior in 
the BI-to-MI transition in one dimension for correlated 
electron systems and establishes a good understanding 
for its ground-state phase diagram. 



VII. SUMMARY 

We have carried out systematic studies using the 
DMRG and TMRG methods to examine the critical be- 
havior of a one-dimensional Hubbard model with an al- 
ternating site potential in the transition from band in- 
sulator to Mott insulator. Based on extensive numerical 
calculations and analytic analysis, we have clarified sev- 
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